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t^J- ' Abstract 

o 

, The computation of long range linear self-interaction forces in string and higher dimen- 

sional brane models requires the evaluation of the gradients of regularised values of divergent 
self-interaction potentials. It is shown that the appropriately regularised gradient in direc- 
tions orthogonal to the brane surface will always be obtainable simply by multiplying the 
regularised potential components by just half the trace of the second fundamental tensor, 
except in the hypermembrane case for which the method fails. Whatever the dimension of 
Oh. the background this result is valid provided the codimension is two (the hyperstring case) 

or more, so it can be used for investigating brane-world scenarios with more than one extra 
space dimension. 



o 
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1 Introduction 

In recent studies of the linear self-interaction of superconducting and other |2|, [| cosmic 
string models in a standard (3+1) dimensional background spacetime, it has been found that 
the relevant divergences can be conveniently dealt with by the application of a simple universal 
gradient formula. This formula can be used to obtain the effective self-force from the appro- 
priately regularized self-field and the precise form of this gradient formula does not depend on 
the physical nature of the fields to which the string may be coupled. These could, for example, 
be of the ordinary (observationally familiar) electromagnetic [|| , and linearised gravitational ||] 
kinds, but might also be of the theoretically predicted dilatonic and axionic kinds |J. 

The purpose of this work is to generalise this gradient formula to higher dimensional branes 
in higher dimensional backgrounds. A particular motivation for doing this is to generalise the 
study of the original kind of brane-world scenario @ not just by dropping the postulate of 
reflection symmetry [||, 01, but by allowing for more than one extra dimension pi, S, 0], s ° as 
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to obtain models consisting of a 3-brane embedded in a spacetime of dimension greater than 5. 
More generally, our results are applicable to any p-brane in a {q + 1) dimensional background 
spacetime, so long as the codimension (q — p) is two (the "hyperstring" case, which includes 
ordinary strings in a (3+1) dimensional spacetime) or more. The method given here does 
however fail (due to long range "infrared" divergences) when the co-dimension is only one (the 
"hypermembrane" case, which includes ordinary membranes in a (3+1) dimensional spacetime), 
but in that case it is not really needed, since the short range "ultraviolet" divergences will have 
a relatively innocuous form that can be dealt with in terms of simple jump discontinuities. 

The reason why one needs a gradient formula is that the relevant forces are typically ob- 
tained as gradients of potentials. Outside the brane worldsheet these potentials will be well 
behaved fields, but in the thin brane limit they will be singular on the worldsheet itself. One 
therefore needs to regularise these fields by some appropriate ultra-violet cut-off procedure, typ- 
ically involving a length scale, e say, that can be interpreted as characterising the underlying 
microstructure. When one has obtained the appropriately regularised potential on the world- 
sheet (which represents the cross-sectional average over an underlying microstructure) the next 
problem is to obtain the corresponding force by evaluating its gradient. The gradient compo- 
nents in directions tangential to the worldsheet can be obtained directly just by differentiating 
the regularised (macroscopic averaged) potential. The trouble is that, since the support of the 
regularised potential is confined to the worldsheet, differentiation in orthogonal directions is not 
directly meaningful. 

To obtain the necessary orthogonal gradient components it is therefore necessary, in principle, 
to go back to the underlying microstructure and perform the ultra-violet limit process again, a 
rather arduous task that has frequently been performed independently by different authors in a 
wide range different physical contexts. 

The present work has been prompted by the observation [|J that, in the case of a string in an 
ordinary 3-dimensional background, the effect of the orthogonal gradient operation will always 
turns out to be equivalent just to multiplication of the regularised value of the relevant potential 
field by exactly half the extrinsic curvature vector K p , which is defined to be the trace 



of the second fundamental tensor of the worldsheet. The recognition of this simple universal 
rule makes it unnecessary to go back and perform the calculation over again every time such 
a problem arises in a different physical context. On the basis of dimensional considerations, 
and of the need to respect local Lorentz invariance, one would expect that multiplication by a 
factor proportional to K p would inevitably be what is required. But what is not so obvious in 
advance is whether the relevant proportionality factor should still always be exactly a half even 
for branes and backgrounds of higher dimensions. 

The present work addresses this question and provides an unreservedly affirmative reply 
whenever < p < q — 1, on the basis of the consideration that whatever the relevant factor 
may be for a generic evolution of a brane of given space dimension p in a back ground of given 
space (as opposed to spacetime) dimension q, this factor must evidently remain the same for 
any static configuration of the brane. 

By restricting our attention to static configurations, for which a rigourous analysis is tech- 
nically much simpler than in the dynamic case, what we have succeeded in showing here is that 
that - provided the codimension {q — p) is two (the "hyperstring" case) or more, so that the 
quantities involved are at worst logarithmically divergent in the "infra red" (which excludes the 
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strongly divergent case of a membrane, but admits the case of a string in 3 dimensions) - the 
appropriate factor is indeed always exactly a half, no matter how high the dimension of the brane 
or the background. However our analysis also shows that this easily memorable result will no 
longer hold in the strongly "infra red" divergent extreme case of a "hyper membrane" (meaning 
a brane with codimension one, so that its supporting worldsheet is a hyper surf ace). What the 
present approach can not do is to provide any evidence at all, one way or the other, about value 
of the factor in question in the opposite extreme case of a zero-brane, meaning a simple point 
particle. 

The article is organised as follows. In Section pi we recall the definitions of the basic geomet- 
rical quantities, in particular the second fundamental tensor, that will appear in our discussion. 
We introduce the Poisson equation, both in the thin brane limit and for regularised configura- 
tions in Section ||| and discuss its solution in Section ||. The cross-sectional average of the field 
and its gradient have a zero-order contribution, in the case of a flat brane, that is evaluated in 
Section [| and we illustrate our formalism on the choice of a canonical profile function in Sec- 
tion §. The cross-sectional average of the gradient of the field vanishes in the flat configuration 
and our goal is to evaluate its value at first order in the curvature. After defining the procedure 
of averaging in Section [7], we study the effect of the bending of the brane on this procedure in 
Section ||. Section || finalises our demonstration by computing the cross-sectional average of the 
gradient of the field. A summary of the notation and some useful integrals are presented in the 
appendix. 



2 Brane embedding geometry 

Since it will play a central role in our analysis, it is worthwhile to recall the definition and some 
basic properties of the second fundamental tensor K^ u p . 

We consider a p-dimensional surface (a p-brane) with internal coordinates a % (with i = 1 . . . p) 
smoothly embedded in a higher g-dimensional background endowed with a metric with 
respect to coordinates (with fj, = l...q). If the location of the surface is given in terms 
of a set of embedding functions x^(a l ), the internal surface metric, also referred to as first 
fundamental tensor, can be expressed as 

dx^ dx v 

^3=9^q-q-- (2) 

Provided that this metric is non-singular (in the sense that the determinant of the component 
matrix does not vanish), as will always be the case for a strictly timelike surface in a Lorentz 
signature background spacetime, or for any surface if the background metric is positive definite, 
then it will have a well defined inverse with components ^ . It can be used to raise internal 
coordinate indices in the same way as the inverse background metric g^ is used for raising 
background coordinate indices. For any such non-singular embedding, j iJ can be mapped into 
the background tensor 

The corresponding mixed tensor 7^ acts on vectors as the natural (rank p) surface tangential 
projector operator, while its (rank q — p) complement 

l»=g$-l$, (4) 
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acts similarly as the corresponding orthogonal projection operator. 

Fields with support confined to the p-surface S^ p \ such as 7^ and can not be directly 
subjected to the unrestricted operation of partial differentiation with respect to the background 
coordinates, but only to differentiation in tangential directions, as performed by the operator 

%=Y^ P , (5) 

where V p is the usual covariant derivative associated with g^ u . The action of this differential 
operator on the first fundamental tensor defines the second fundamental tensor of the p-surface, 



according to the specification |1C] 



AV = i) w V/ • (6) 

As a non-trivial integrability condition, it satisfies the generalised Weingarten symmetry condi- 
tion 

AV = A,/ . (7) 
It is also p-surface orthogonal on its last index and tangential on the other two, that is, 

K^rf a = =_L* KJ> , (8) 

so that it has only a single non-identically vanishing self-contraction, namely the extrinsic cur- 
vature vector introduced in equation ([I]). 

3 Regularisation of brane supported source distribution 

Since we are essentially concerned with the treatment of "ultraviolet" regularisation, the effects 
of long range background curvature will be unimportant and can thus be neglected. As the 
framework for our analysis, it will be sufficient to consider a q-dimensional background space 
with a flat, static, positive definite metric g^. It is therefore possible to choose a linear (but 
not necessarily orthogonal) system of space coordinates so that the background metric will 
be constant, that is, dg^/dx 9 = 0. 

This choice has different implications. First, the tangentially projected covariant derivative 
^ will take the simple form 

V -Y— -Q ^ — — (9) 

dxP da i da r W 

Second, in the linear system (for, example Maxwellian or linearised gravitational), with which we 
are concerned here, each tensorial component decouples and evolves independently like a simple 
scalar field. Thus, it will be sufficient for our purpose to concentrate on the prototype problem 
of finding a scalar field, let us say 4>{x} (introducing the systematic use of curly brackets to 
indicate functional dependence), that satisfies the corresponding generalised Poisson equation 
to which the linear wave equation will be reduced in the static case to be dealt with here. Third, 
the background spacetime is assumed to be static so that the generalised Laplacian V M V M is to 
be regarded as as the static projection of what would be a Dalembertian type wave operator in 
a more complete (p + l)-dimensional spacetime description. 
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3.1 Generalised Poisson equation and its solutions in the thin brane limit 

In unrationalised units for, let us say, a repulsive scalar source distribution p{x] the generalised 
Poisson equation takes the form 

v^<f> = -n l9 - 1] P {x}, (io) 



where the Laplacian operator is simply given by 

dx^dx u ' 

(<?-!> 

latter is given by the well known formula 



where iV 9 11 is the surface area of a unit (q — l)-sphere in the g-dimensional background. The 



n"-" = ^r, (12) 

r{«/2}' v ' 

in which T{z} is the usual Eulerian Gamma function (as specified in the appendix) which 
satisfies the recursion relation T{n + 1} = nT{n} and has the particular values r{l} = 1 and 
T{l/2} = yfR. Thus in particular we shall have ftM = 2, fiW = 2vr, fiPl = 4vr, ft! 3 ! = 2vr 2 and 
nW = 8vr 2 /3. 

The motivation for working with a p-brane model of a physical system is of course to be able 
to use an economical description in which as many as possible of the fields involved are specified 
just on the supporting surface. Evidently, this requires less information than working with fields 
specified over the higher dimensional background. However, when long range interactions are 
involved then, as a price for such an economy, divergences are to be expected if the source 
fields are considered to be strictly confined to the supporting p-surface. To be more specific, 
the brane supported source distribution p{cr} (using a bar as a reminder to indicate quantities 
that are undefined off the relevant supporting surface) will correspond to the background source 
distribution 

p{ x } = J 6 [q] {x,x{a}}p{a}dS lp] , (13) 

where <5 M {x, x{a}} is the g-dimensional bi-scalar Dirac distribution. This is characterised by 
the condition that it vanishes wherever the evaluation points x M and x^{a} are different while 
nevertheless satisfying the unit normalisation condition 



5 {x,x}dS =1. (14) 

We use the notations d<S W and d5 M respectively for the p-dimensional brane surface measure 
and g-dimensional background space measure. The latter is given by an expression of the form 

dS W = M 1/2 d M z, (15) 

where the g-dimensional volume measure factor 

\g\ 1/2 is 

just the square root of the determinant 
of the (positive definite) matrix of metric components g^ u . Similarly, we have 

d<S [pl = | 7 | 1/2 d W a, (16) 
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where the p-dimensional volume measure factor I7I 1 / 2 is the square root of the determinant of the 
matrix of surface metric components, as given by (|2|). The normalisation condition (|l4j) is such 
that the measure factor means that the Dirac distributions specified in this way will transform 
as an ordinary bi-scalar, a property that distinguishes it from the corresponding Dirac type that 
would be given by the commonly used alternative convention in which the measure factor is 
omitted. This specifies distributions that behave not as ordinary bi-scalars but as weight one 
bi-scalar densities. Of course it does not matter which of these conventions is used if, as will be 
done below, the coordinates are restricted to be of the standard orthonormal type with respect 
to which the metric just has the Cartesian unit diagonal form. 

For the source distribution given in the thin brane limit by ([ilf), the corresponding solution 
of the linear field equation ( |l0| ) is given in terms of a Green function by the expression 

4>{x} = J G l " ] {x, x{a}} p{a} dS M , (17) 

while the corresponding scalar field gradient - which is what determines the resulting force 
density - will be given by 

r dG lq] M 
M= / -^{x,x{a}}p{a}dS [p] . (18) 



dxv J dx^ 

The Green function G M {x, x{a}} satisfies the Laplace equation 

V^{i,iW} = -n lq ~ 1] 5 l " ] {x,x{a}} (19) 

subject to the appropriate boundary conditions, which in the standard case for an isolated system 
will simply consist of the requirement that the field should tend to zero at large distances from the 
source. Like the Dirac "function" , it is a bi-scalar (not a bi-scalar density) that is well behaved 
as a distribution, but singular as a function wherever the two points on which it depends are 
coincident. Note that the Green function is determined up to an harmonic function, that is a 
solution of the homogeneous Poisson equation V^V^G 19 ' {x, x{a}} = 0, with the appropriate 
boundary conditions. In particular G^ is determined up to an additive constant that induces a 
global shift in the potential and is not relevant for the computation of its gradient. 

3.2 From the thin brane to a regular thick brane 

Although convenient for many formal calculational purposes, the drawback to the use of the 
strict thin brane limit as represented in terms of such singular Dirac and Green distributions 
is that the ensuing value of the corresponding field cj) and its gradient dcp/dx 11 will also be 
singular on the brane p-surface, which unfortunately is just where one needs to evaluate them 
for purposes such as computing the corresponding energy and force density. 

The obvious way to get round this difficulty - without paying the cost of introducing extra 



degrees of freedom - is to replace (13) by an expression of the form 



~ P{e) {x} = J vf e} {x,x{a}}p{a}dS lp] , (20) 

in which the strict Dirac distribution # M {x, x} is replaced by a regular bi-scalar profile function, 
P { } {x, x} say, characterised by some sufficiently small but finite smoothing length scale, e say, 
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which is subject to the normalisation condition 

/p;: ] } {x,x}d5 M =l. (21) 

In order to respect the geometric isotropy and homogeneity of the background space, this regular 
"profile function" X>j* {x,x} can depend only on the relative distance, s say, that is defined - 
in the flat background under consideration - by 

s 2 = 9^ - x»){x v - x v ) , (22) 

according to some ansatz of the form 

vf t} {x,x} = V 1 ^}. (23) 

The functionally singular Dirac distribution can be considered as the limit of such a function 
as the smoothing lengthscale tends to zero. Such a description corresponds to what may be de- 
scribed as "fuzzy" brane model characterised by an effective thickness whose order of magnitude 
will be a function of the value of the regularisation parameter e. For some kinds of physical 
application one might wish to adjust the regularisation ansatz so to make it agree as well as 
possible with the actual internal structure (if known) of the extended physical system under 
consideration. The spirit of the present work is, however, rather to treat the regularisation 
parameter just as a provisional freely variable quantity for use at an intermediate stage of an 
analysis process whereby it would ultimately be allowed to tend to zero so as to provide a strict 
"thin" brane limit description. In order for our results to be meaningful for this latter purpose, 
it is important that they should be insensitive to the particular choice of the smoothing ansatz 
that is postulated. We shall take care to verify this 



4 Field solution for regularised source distribution 

For any source distribution of the regularised form (20), the function, C/|^{s 2 } that is the solution 
(with the appropriate boundary conditions) of the corresponding generalisation 

v^;: ] } { S 2 }=-^ [9 " 1] < , } {^ 2 }, (24) 



of fll9| ) will act as a generalised Green function. The solution of the Poisson equation ( |i0| ) for 
the generic smoothed source distribution ( p0|) will thus be given, using the natural abbreviation 



s{x, a} = s{x, x{a}} , (25) 

by 

to = / C W*. M <±s lP] • (26) 

The gradient of this scalar field is evidently given by 

= / S?^' a}} ^ ' (27) 
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in which the integrand can be evaluated using the formula 

as'", ds'", 

where = x^ — x^. One can evidently go on in the same way to obtain the second derivative 
in which the integrand can be evaluated using the formula 

d 2 g\ q \ dg l ;\ d 2 g [ * 

'" -^-TT + ^^ujt^- (30) 



dx»dx» ~ yflu ds 2 ^ " "(ds 2 ) 2 ■ 



It is to be noted that, by contraction of (|30|), the Laplace equation (|4|) for takes the form 



,! ; A>i --^«- 2 p;: ] } {, 2 }. (3i) 



ds 2 I ds 2 



If one introduces the monotonically increasing function zj^{s 2 } by 



s 



T;\{s*} = tt H ' ui- L V"'{u 2 }du, (32) 

U} Jq {e} 

then it will satisfy 



ds 2 2 

By substituting into (31), it follows that we shall have 



dT M r> lq ~ 1] 



ds 2 2s<? ' 

and hence 

{«} _ O} L J {<=} L J 

(ds 2 ) 2 ~ 4si+ 2 4s 2 



(34) 



(35) 



The normalisation condition ( |2~l| ) implies by ( |32"| ) that the function T { ^ } {s 2 } should satisfy the 
condition 

<:> 2 }-i (36) 

in the limit s 2 — > oo, for all e. It follows that any admissible ansatz for the smoothing profile 
will entail an asymptotic behaviour characterised by 

d *, x 1 

— ~ (37) 

ds 2 2s"? v ; 
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as s 2 — ► oo. The corresponding asymptotic formula for the regularised Green function Q^{s 2 } 
will therefore be given, as long as q > 2, by 

as s —* oo. In the particular case (7 = 2, the integral over s 2 is marginally convergent, with 
asymptotic behaviour characterised by 

<> 2 } ~ ~\Hs 2 } (39) 

as s 2 — > oo. 



5 Evaluation for flat brane configurations 

Before considering the effect of the curvature of the brane, it will be instructive to apply the 
preceding formulae to the case of a flat p-brane configuration supporting a uniform source 
distribution. 

The uniformity condition means that the value of p is independent of the internal coordinates 
a 1 , and the flatness condition means that the embedding mapping a 1 i— > ^{a 1 } of the brane p- 
surface will be expressible in a system of suitably aligned Cartesian background space coordinates 
{x^} = {z l ,f a } (with fj, = 1 . . . q, i = 1 . . . p, a = p+ 1 ... q) in terms of a corresponding naturally 
induced system of internal coordinates a % simply by 

2* = a\ r a = 0. (40) 

With such a coordinate system, the background metric components will be given by 

9ij = liji 9ia = 0, g a b =-L Q 6 , (41) 

where jij and _L a b are unit matrices of dimension p and (q — p) respectively. 

The intrinsic uniformity of the configuration ensures that physical quantities depend only 
on the external coordinates r a , and the isotropy of the system ensures that scalar quantities can 
depend only on the radial distance r from the brane p-surface given by 

r 2 =± ab r a r b . (42) 

In particular, the distance s from a point with coordinates on the transverse {q — p)-plane 
through the origin (that is, with a 1 = 0) to a generic point with coordinates a on the flat p-brane 
locus will be given simply by 

where 



s 2 {x,a} = r 2 + a 2 , (43) 



a z = 7y aV . (44) 

It can be seen that for a given fixed value of the source density p on the brane, the prescription 
(|20| ) will provide a radial dependence of the form 

p{r 2 }= ± V l ^{r 2 }p, (45) 



9 



where the dimensionally reduced radial profile function {f 2 } is given in terms of the 

original g-dimensional profile function T>^{s 2 } by the expression 

V 9 f {r 2 } = d p - 1] [°° a^V^la 2 + r 2 } da . (46) 

This integral will always be convergent as a consequence of the rapid asymptotic fall-off condition 
that must be satisfied by the function P M {s 2 } in order for it to obey the normalisation condition 
(21), which implies that 

Q lq - p - 1] J™ ri-v- 1 X>V 2 } dr = 1 . (47) 

Since (j37|) implies that the associated regularised Green function C/j^{,s 2 } will satisfy the 
comparatively weak asymptotic fall off condition (|38|), the field (f), , will be given by 



^ } {r 2 }=- p X?{r 2 }, (48) 
using (p6[), where the dimensionally reduced radial Green function X ^|' } P '{^ 2 } is defined by 

X? {r 2 ) - ^ jf o*- 1 ^ {- 2 + da , (49) 

in the limit where A —* oo, a limit that will be convergent only if g > p + 2. 

In the prototype case 0, || of an ordinary string with p = 1 in a background of space 
dimension (? = 3, the latter is only marginally too low compared with the dimension of the brane 
surface to achieve convergence. In that case, as in the more general hyperstring case, that is to 
say whenever q = p + 2, it is still possible to obtain a useful asymptotic formula by taking the 



upper limit A to be a finite "infra-red" cut off length. It can be seen from (38) that in this 



marginally divergent codimension-2 case the resulting asymptotic behaviour will be given by 

n lq ~ 3] . r a 



»« ^~—^{wttA< (50> 

as A — > oo. 

One might go on to try obtain an analogous formula for what, as far as "infra-red" divergence 
behaviour is concerned, is the worst possible case of all, namely that of a hypermembrane (or 
"wall") meaning the hyper-surface supported case characterised by q = p + 1. However in this 
codimension-1 case, the divergence will have the linear form 

which is particularly sensitive to the choice of cut-off A and hence is not very useful. 

For physical purposes the potential </> {e} will typically be less important than its (not so 
highly gauge dependent) gradient, which is what will determine the relevant force density. In 
the flat and uniform configuration considered in this section, the gradient d(j)^ / dx^ is not liable 
to the infrared divergence problems that beset the undifferentiated field when the background 
dimension p is not sufficiently high compared with the brane dimension q. It is evident from the 
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uniformity of the configuration that the gradient components in directions parallel to the brane 
p-surface will vanish, that is, we shall have d<j),,/dz l = 0, while starting from the expression 
(p7|), it can be concluded that the orthogonal gradient components will be given by 



i^{r 2 } = P^{r 2 } (52) 

where, by Eq. (j34|), even in the cases for which the undifferentiated integral (fls|) would be 
divergent, the required derivative is given by a well defined formula of the form 

^{^ = -r a ^I;V}, (53) 
in which the function J^'^ {r 2 } is defined as the integral 

, „ roo l lq] {r 2 +a 2 } 

which will always converge by virtue o f (|36|) . 

In a similar manner, starting from (|29|), it can be seen that the corresponding second deriva- 
tive components will be given in terms of the same convergent integral J7i*' Pl {r 2 } by 

d% } 2 _ d 2 x g;*f 2 

Q r aQ r b "L 7 * } ~~ P Q r aQ r b { T J (55) 

in which, by the definition (|49|), it can be seen that we shall have 

d 2 ± g lq ' p] , „ /-a / dg l9] d 2 g l9] \ 

■rff > - M i ^ ■# + 2r -id#) f + ^ «° ■ ^ 

in which it is to be recalled that that the use of curly brackets (as in the expression {a 2 + r 2 } 
at the end) indicates functional dependence (not simple multiplication). Using ( |35| ) to express 
the second derivative in the integrand as 

(ds 2 ) 2 2r 2 ds 2 4r 2 ' r 2 da 2 I ds 2 J ' W j 

and observing that the asymptotic behaviour (|37| ) ensures that the last term in ( |57| ) will con- 
tribute nothing to the integral on the right of (|56|), we finally obtain an expression of the 
convenient form 



8 2 ± G l "' pl 



^^{V} = -^<r ^-p)^-^ 6 j47 ] {r 2 } ; (58) 

in which J, Z?|*'*' {r 2 } is the dimensionally reduced profile function given by the integral (|4(f). 
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6 Illustration using a canonical profile function 



Our main results will be derived for any generic smoothing ansatz, without restriction to any 
particular form for the profile function V^{s 2 }. However as a concrete illustration, it is in- 
structive to consider the case of what we shall refer to as the "canonical" regularisation ansatz 
whereby the g-dimensional Dirac delta "function" and the associated Green function 

G M (x, x) are considered as the respective limits - as the regularisation parameter e — > - of 
one-parameter families of smooth bi-scalar functions 

^ ] } {*,*} = ^ } { S 2 } , and gf e) {x,x} = G^s 2 } (59) 

that are characterised in terms of the squared distance s 2 {x, x} and of the smoothing lengthscale 
e by the prescription 

5 u } i s } = ^r( s + e ) ■ ( 6 °) 

This canonical profile function evidently satisfies ^^{s 2 } — > when e 2 — > oo for s 2 ^ 0, 
as well as the unit normalisation condition (21) using the fact that dSM = O 1 " s'^ds. By 



expressing the canonical profile function as 

q/2 



*: } {s 2 } 



2 2- q d 



-s 



O 19 - 11 ds 2 



,s 2 



s 2 + e 2 



it can be seen that the resulting integral (|32"|) will simply be given by 



(61) 



<,f^(^)' /2 - («) 

It is then straightforward, by integrating (|34|), to show that the corresponding regularised Green 
function, satisfying by the Poisson equation (p4|), will be given by the canonical Green function 
defined by 

Gl } {s 2 } = j^(s 2 + e 2 )-^\ (63) 

as long as q > 2. As explained above, it is determined up to a constant that is fixed by the fall 
off requirement at infinity which is mandatory for the convergence of the integrals considered 
below but that will not be relevant for the gradients to be computed later on. In the marginally 
convergent case where q = 2, one finds that 



For a specific physical system this particular choice is not necessarily the ansatz that will 
be most accurate for the purpose of providing a realistic representation of its actual internal 
microstructure. However this "canonical" prescription - as characterised by the formulae ( |6"2| ) 
and ((6^) - is distinguished from other conceivable alternatives by the very convenient property 
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of having an analytic form that is preserved by the dimensional reduction procedure. The 
associated radial profile function is simply given by 

V [q ' V Ur 2 \ - - Q / dv (6*)) 

U U} ^ ' n h-H ( e 2 + r 2)( g - P +2)/2 7 (1 + n 2)( g +2)/2 QU " 

Since we shall always have g > p — 2, the latter integral will always be convergent. Using the well 
known properties (see (118) in the appendix) of the Euler Beta function B{p/2, (q — p)/2 + 1} 
it can be seen that in this canonical case the radial profile function will be given simply by 

-Lptw | ,.2 i . A -h j ;2 i 

It can similarly be seen that, since we always have q > p, the integral (|54j ) will also be convergent, 
and (again using the formulae in the appendix) that the result will be expressible in terms of 
the Beta function as 

^> {r} = 2(^+^)to- P )/^ \2'^r (67) 

According to ( fl9| ) the corresponding dimensionally reduced Green functions will have a 



canonical form 

^Vl^^Vl (68) 
that will similarly be given whenever they are well defined - that is whenever q > p + 2 - by 
expressions of the form 

n b-i] , 



e ^ = ^-2)^ + ^^ B \i-"-r- 1 i' (69> 

so that, in terms of the canonical Green function ([63]), they take the form 

V ? f {r 2 } = ^-G l r x P] {r 2 } . (70) 

{e} L J QW-P-1] {e} L J V ' 

The form of the Green function, unlike that of the canonical profile function, is thus not exactly 
preserved by the dimensional reduction process but the only difference is an overall volume 
factor. The preceding expression is only valid when q > p + 2, but in the marginally convergent 
hyperstring case where q = p+2 one can use the expression (|64| ) for the canonical Green function 
to obtain the asymptotic relation 

n [q - 3] . r a 



C i r } ~ 7^ — T ln -^== ( 71 ) 

W 1 ; ( g -2) \V^T?J 

as A — > oo. It is to be noted that the logarithmic term in this formula and in ( |50| ) will be 
expressible by d64|) in the form 



ln 



{^M> 2 >+>»{!} < 72 > 



in which the first term is independent of the cutoff A. This means that in the asymptotic limit 
the dependence on r will drop out, so that we shall be left simply with 

,[9-3] 



-L„[8.«-2] r 2t ^ 



as A —* oo. 
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7 Source weighted averages in flat canonical case 



For the purposes of a macroscopic description, what really matters is not the detailed field 
distribution in a smoothed microscopic treatment as developed in the previous sections, but 
only its effective averaged values. For instance, in the example of the scalar field </> {£} {r}, we 
shall be interested in its cross-sectional average defined as 

(^ {e} ) = n [ "- p - 1] f™r*-*-% } {r}w{r}dr, (74) 

where w{r} is a weighting factor dependent only on the radial coordinate r subject to the 
normalization condition 

d q ~ v ~ 1] r r q -*- x w{r} dr = 1 . (75) 



o 

Let us first consider the flat uniform configurations studied in Section [5|. Whatever the 
explicit form of the function w{r}, it is obvious that the isotropy of the configuration will 
ensure that the average of the field gradient ( |52[ ) will cancel out, meaning that we shall have 
d(f)^/dr a ^ = 0. Using the general relation 

™ X 1 ± ab (4r), (76) 



f{r}/ q-p \f{r} 

which applies to any purely radial function f{r}, it can thus be seen from (|55| ) and (|5q) , using 
the expression (||) for the radial source density distribution, that 

(d^p) = -J—£( p w) ±ab > (77) 

in which the value of (p^y^ a l so depends on the choice of the weighting function. It is to be 
remarked that (as a useful check on the preceding algebra) the trace this equation gives back 
the original Poisson equation (|56|). 

The kinds of energy and force contributions that are relevant in physical applications will 
typically be proportional to the product of the linear field <p with the corresponding source term. 
This implies that the weighting of the field </> {e} considered here should be proportional to the 
source term p {e} . For a flat configuration, it can thus be deduced from (^) that the appropriate 
weighting factor should be explicitly given by the formula 

-M=X>Vl, (78) 
and hence that the average of the source term will be given by an expression of the form 



= pd q - p - 1] I X;> 2 } " r q - p ~ 1 dr , (79) 
Jo 1 {£} J 

which will be valid, whatever the profile function, as long as the configuration under considera- 
tion is flat. The corresponding expression for the average of the field <p will have the form 

<0 {E} ) = pV lq - p - 1] f { r ^ {r 2 }r q ~ P ~ l ^ . (80) 
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In the special canonical case introduced in Section ||, the profile factor will be given by 
|), so that the (always convergent) integral (79 will be obtainable explicitly - as shown in the 
appendix - in terms of the Beta function. This leads to an expression of the form 

P< >) = ^ + 2 1 J- (81) 

' 2f2 I 2 2 J e q ~P 

for the average source density. With further use of the formulae in the appendix, this can be 
written even more explicitly as 

\ 1 (q-p) 2 (q-p + 2)[T{(q-p)/2}} 3 p 

{E>/ lOvr 2 ^ (q-p+l) T{q-p} e*-P ' 



The simplest application of this result is to the case of a hypermembrane or "wall", meaning 
a hypersurface forming brane as characterised by q — p = 1, for which the formula ( |82[ ) just 
gives ^/°{ e }^) = 37r / o/(32e). For the case q — p = 2, which includes that of a string in an ordinary 

3-dimensional space background, the corresponding result is (p{ e y) = p/(37re 2 ). 

In terms of the canonically reduced Green functions ( ]68| ) , it it can be seen that the average 
of the the regularised field (^ {e} will be given in the flat case by a prescription of the form 

(* w ) = pd q ' P ~ 1] r ^{r^^r^dr. (83) 



Whenever the expression ( |70|) is valid, that is whenever q > p + 2, the formulae in the 
appendix will provide a result expressible as 

1 {q ~ p) --' "(V^)^- m 



(85) 



{E} / 2 (q - p - 2) n lq ~ p ~ 1] I 2 ' 2 J ei~P~ 2 ' 
which can be rewritten in terms of V functions as 

* p/2 (q-p) [T{(q-p)/2}f p 



{A 



2 {q-p-2)T{q/2}T{q-p} ei-v- 2 ' 



The preceding formula will no longer be valid, but an analogous though weakly cutoff de- 
pendent expression will still be obtainable in in the marginal hyperstring case, when q = p + 2. 
In this case, instead of (70), the effective dimensionally reduced Green function will be given 
by the simple asymptotic formula (73) in the large value limit for the relevant infrared cutoff 
A. The ensuing asymptotic formula for the averaged value of (r) in a hyperstring is thus 
immediately seen to be given by 



pri q - 3] (A 

i'l/ ~2{,/ 2) I ( 



(86) 



as A — > oo. 
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8 Allowance for curvature in the averaging process 



We are now ready to tackle the main problem that motivates this study, which is the evaluation 
of the effect of curvature of the brane p-surface. 

In order to do this, instead of the exactly flat configuration characterised by the coordinate 
system (40) and the particular form of the metric components ([Il|), we need to consider a 



generically curved p-surface adjusted so as to be tangent to the reference flat p-surface at the 
coordinate origin. 

The validity of the "thin brane" limit description with which we are concerned depends 
on the requirement that the relevant macroscopic lengthscale A (that was used to specify an 
appropriate "infra-red" cut off, in cases for which the convergence condition q > p + 2 does not 
hold) should be large compared with the microscopic brane thickness scale e. Subject to the 
requirement that 

a = ^ « 1 (87) 

is satisfied (either because A is very large or because e is very small), it will be sufficient for 
our present purpose to work at linear order in a. This means that we need only consider up to 
quadratic order in a 1 for an expansion where the brane locus will be specified in the relevant 
neighbourhood - of dimension large compared with e but small compared with A - by 



i^- W [1 + o{a)}\ , 



for an arbitrary set of constant coefficients Kij a . In this expansion, we have used the standard 
meaning of o{5} denoting a quantity negligible with respect to 5, that is satisfying o{5}/5 — ► 
in the limit 5 — > 0. It can be easily verified that, as our notation (|8g| ) suggests, the constants of 
the development are exactly given by the non vanishing components of the second fundamental 
tensor K^ v p ', defined in @ at the origin, with respect to the orthonormal background coordinates 

In the simplest cases, the condition ( |S7|) might be satisfied simply by taking A to be arbi- 
trarily large. In specific applications an upper bound on the admissible magnitude for A will 
commonly be imposed by the physical nature of the environment in which upper limits might 
be provided by lengthscales such as those characterising the gradients of relevant background 
fields or those characterising the separation between neighbouring branes. However, an upper 
limit that must always be respected will be provided, except in the strictly flat case, by the cur- 
vature lengthscale defined by the inverse of the magnitude of the mean curvature vector whose 
components are given by the formula 

K a = K 4 ia . 



As our notation suggests, these components K a are identifiable with respect to the orthonormal 
background coordinates ( |40| ) as the non-vanishing components of the complete curvature vector 
given by (||). In order to ensure that none of the coefficients Kij a exceeds the order of 
magnitude of the inverse lengthscale A -1 , the latter must thus be subject to the limitation 

A~ 2 > K a K a . (90) 

In particular, our formalism will break down near a kink or a cusp of the brane worldsheet. 

The curvature of the worldsheet makes the evaluation of the physically appropriate weighting 
factor for the definition of effective sectional averages across the brane rather more delicate than 
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in the flat case dealt with in the preceding section. The kind of average that is ultimately relevant 
is a source density weighted mean taken not just over a (q — p)-dimensional orthogonal section 
through the brane but rather over the q-dimensional volume between neighbouring {q — p)- 
dimensional orthogonal sections through the brane. This distinction is irrelevant for a flat brane 
for which neighbouring sections are exactly parallel. However for a slightly curved brane, there 
will be a relative tilt between nearby orthogonal cross sections, which means that the effective 
weighting factor w{r} appropriate for the purpose of averaging over a (q — p)-dimensional cross- 
section will no longer be exactly proportional to the density p{r} as in the flat case. This needs 
to be adjusted by a correction factor allowing for the relative compression or expansion of the 
relevant (/-volume due to the relative tilting. By summing over the different directions in which 
the tilting can occur, it can be seen that the weighting factor will be given, at first order, by 



Pi x{ r } 

w { r } = lillli. [i _ Ka r a + o{a}] . 
P 



(91) 



To obtain the corresponding generalisation of the formula (78) for the flat case, we must 
allow for the fact that the original formula ( |45| ) for the source density distribution also needs to 
be corrected to take into account the effect of the curvature. Using (^) to expand the distance 
function ( p2| ) as 

S 2{ X) a } = r 2 + a 2 - ± ab K tJ a r b a 1 ^ [1 + o{a}} , (92) 
we obtain the Taylor expansion 



dV lf {r 2 + a 2 } 



^ } {* 2 }=^ } {^+- 2 } 



± ab K^a 1 ^ [1 + o{a}} 



da 2 



(93) 



for the profile function fl23|). Since we are neglecting corrections of quadratic and higher order 
in the curvature, it can be seen that when the integration is carried out in two stages of which 
the first is just to take the integration over the p-sphere of constant radius a centered on the 



origin, it is possible to replace the surface element dS [p] in ( |20| ) by fi 1 " iJ cr p-1 dcr as in fl46|). 
Furthermore, in the anisotropic contribution proportional to a % a 3 it will be possible to replace 
crV^dc? 1 by the spherically averaged equivalent j^^ P 11 a p+1 da /p. Therefore, to allow for the 
effect of the curvature at first order, the density distribution (|45|) needs to be replaced by 



Jp-i] 



P U} {r} = P 



V 



li,p] 



{<=} 1 s 



[K a r a + o{a}] / a^da 2 

V Jn da z 



(94) 



in which the second term can be evaluated using an integration by parts, using the definition (|4 
for ± J)^'^ {r 2 } and the final result takes the simple form 



P {c} {r}=P V {A {r } 



1 + l -K a r a + o{a} 



(95) 



Thus, it can be seen that the geometric curvature adjustment factor in (91) is partially canceled 
by the density adjustment factor in (95) to give the corresponding curvature adjusted weighting 
factor 

1 



w{r\ = V {e} {r } 



1 - -K a r a + o{a} 



(96) 
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We note that as a corollary of this partial cancellation, the first order effect of the curvature 
is entirely canceled in the correspondingly weighted density function, which will be given simply 
by 

^{r}p {e} {r} = p[ ± V [ ^{r 2 }} 2 [l + o{a}} . (97) 



Thus to first order in the curvature, the average density will still be given by the same form (|7 
as in the flat case. 



9 Allowance for curvature in the potential gradient 



In the flat configuration, the average of the field gradient (^d(j)^ /dr a j , vanishes. Thus, un- 
and (d 2 4>, } /dr a dr b \ for which the dominant contributions, respectively given by 



like 



Eqs. (74) and (|77|), are of zeroth order in the curvature, it will have a dominant contribution 
that is of linear order in the curvature. 

Starting from the general expression (26) for the field </> {e} in terms of the generalised Green 
function, the strict analogue of the integral (|9^) for the source density /? {e} is given, after Taylor 
expanding the Green function, by 



\ e} {r} = pQ 



b-i] 



a 



a 2 d ^ } {r 2 + a 2 } 
p da 2 



K a r a 



[l + o{a}\. (98) 



The truncation of the integration at the finite infra-red cut-off A is only needed for the treatment 
of the marginally divergent case for which q = p + 2, but makes no difference, at first order in 
e/A. In the convergent cases q > p + 2, the result will be effectively the same as would be 
obtained simply by taking the limit A — > oo. The corresponding expression for the gradient is 
found, after an integration by parts, to be 



{<=} 



dr a 



pQ 



[p-i] 



a- 



p - l da 



(2 + K b r b ) r a - jK, 



2 + a 2 } 

[1 + o{a}} . (99) 



The appropriately weighted potential function required for the evaluation of the relevant 
average can now be seen from ( |96| ) and (|9l| ) to be given, after another integration by parts, by 



.[9-1] -ly-,[8.P] 



™{r} 4> U} M = fin—' X>> 2 } / ^tog^ii* + [1 + o{a}} 



-pSl 



[9-1] 



,Af 



(100) 



When q > p + 2, the asymptotic behaviour fl38]) implies that the boundary contribution tends 
to zero, while the integral converges toward the radial Green function, as defined in (|4~9| ) so that 
we obtain an expression of the simple cut-off independent form 



w{r} 4> {c} {r} = p ^V [q { ^{r 2 } ^^{r 2 } 



In the marginal case q = p + 2, we obtain an expression of the asymptotic form 
w{r} <P {t} {r} = p ^V l ™ *'{r 2 } 



In 



K n r a 



2( g -2)2 



(101) 



(102) 
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in which the boundary contribution at the end will provide a finite curvature adjustment term. 
However, since this adjustment term is an odd function of the radius vector r a it will still provide 
no net contribution to the integrated average, which will thus be given simply by 



{<=} 



pn^'P^Mr^ir 2 } ^{r 2 } [1 + {a}] , (103) 



which to first order in a is just the same as in the corresponding formula ( p0[ ) for the flat case. 

When we go on to work out the analogously weighted potential gradient distribution, it can 
be seen that the term with quadratic radial dependence in (|99[) will cancel out, leaving 



r -i -n lp - 1] Vi^r 2i f p-lj 

w\r\ — — = nil D,, r cr^ da 





a 2 

2r a K a 


/o 


V 



ag [q \{r 2 + a 2 } 

— ^ [l + o{a}]. (104) 



Since the first term in the previous integral is an odd function of the radius vector r a , it provides 
no net contribution to the corresponding integrated average, which will therefore be given simply 
by 

roo 



Q r a 



>d"- p - 1] K a nr^&r^{r*}H [q > v \r*}{l + o{a}] , (105) 
Jo ic> 



n [q - p] {r 2 } = / a^da U} \ = . (106) 



where 

p jo do -2 
Using an integration by parts, this integral can be rewritten in terms of the dimensionally 
reduced Green function G^y^'^ir 2 } defined by ( [49]) in the form 

[q,p] , 2l 1 _L [«,p] , 2l ^ ^ P r 2 , A 2l / 1n «\ 

W i r } = 2 G {e} V } ^^^} {r+A} - (107) 

As in ( |1Q0|) , the boundary contribution at the end will vanish in the large A limit in consequence 
of the asymptotic behaviour ( |3"8| ) provided the convergence condition q > p + 2 is satisfied. In 
the marginally divergent hyperstring case q = p + 2, the boundary term will tend to a limit 
having a finite value (namely —l/2p 2 ) that will still be negligible compared with the first term 
in ( |107| ). Thus, whenever q — p > 2 we shall always have a relation of the form 

n [qM {r 2 }= l -G [ ^{r 2 }[l + o{a}]. (108) 

The only case to which this formula does not apply is that of a hypermembrane, meaning 
the hypersurface supported case characterised by q = p — 1, for which the boundary term will 

We can still derive an asymptotic relation 



be linearly divergent, like the corresponding the integral for P ' {r } as characterised by (|5l|). 



n lq - q - 1] {r 2 } ~ ^yXr ^ ( 109 ) 

as A — > oo even in this extreme case, but its utility is limited by the strongly cut off dependence 
of the quantities involved. It can be seen that it approaches agreement with the generic formula 
(|108j) when the space dimension q is very large. 
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When (|108|) is substituted back into ( |l05l) one obtains a result that can immediately be seen 
by comparison with the formula ( |103| ) for the averaged potential ^0{ e }^) to be expressible in 
terms of the latter by 

1°'^' \ l K a U^)[l + o{a}] , (110) 



\ dr a 2 a \ {e} 

a simple and easily memorable result whose derivation is the main purpose of this work. 

10 Discussion 

The preceding work shows not only that the relation ( |110| ) holds as an ordinary numerical 
equality for the strictly convergent cases for which the co-dimension, (q — p), is greater than 2, 
but also that it holds as an asymptotic relation for the marginally divergent case of a hyperstring 
with q — p = 2, including the previously studied examples |l], ||, || of application to a string with 
p = 1 in the ordinary case of a background with space dimension q = 3. 



It is also clear from the preceding work that the formula ( 110 ) is not applicable to the 



case of hypermembrane, for which the co-dimension is given by q — p = 1. In this case, it 



can be seen from ( 109 ) that the factor 1/2 in ( |110D should in principle be replaced by a factor 
(p— l)/2p. However it is debatable whether any useful information can be extracted in this case 
since the result is strongly dependent on the cut-off due to the "infra-red" divergence, and in 
practice this issue is of little importance because the particularly good ultraviolet behaviour of 
the hypermembrane case makes it easily amenable to other, more traditional, methods such as 
the use [|ll]] of Israel Darmois type jump conditions. 

Dropping the explicit reminder that higher order adjustments would be needed if one wanted 
accuracy beyond first order in the ratio of brane thickness to curvature radius, the relation ( |110| ) 
translates into fully covariant notation as 

^V^ {£} ) = V(^ {£} ). (Ill) 



It is to be emphasised that the generic relation (111) is a robust result, whose validity, whenever 



the co-dimension is 2 or more, does not depend on the particular choice of the canonical regular- 
isation ansatz on which the explicit formula (85) for \</> {£ } ) was based. It also shows that ( |111| ) 



is valid for any alternative non-canonical regularisation ansatz, whose mathematical properties 
would be less convenient for the purpose of explicit evaluation, but that might provide a more 
exact representation of the internal structure for particular physical applications in cases where 
information about this actual internal structure is available. 

The question that remains to be discussed is how this result generalises from static Poisson 
configurations, as considered here, to dynamic configurations in a Lorentz covariant treatment 
for which the Laplacian operator would be replaced by a Dalembertian operator. 

A priori, considerations just of Lorentz invariance imply that the only admissible alternative 



to the formula ( 111 ) would be a formula differing just by replacement of the factor 1/2 by some 
other numerical pre-factor. However, it must be taken into account that this numerical coefficient 
must match the coefficient that applies to any static configuration. It can thus be deduced from 
the postulate of consistency with what has been derived here that the formula (|111| ) should 
be applicable to any p-brane with p + 1-dimensional timelike worldsheet in Lorentzian q + 1- 
dimensional background spacetime, for all values of p and q for which our present derivation 
applies, that is to say for q — 2 > p > 1. 
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There are two interesting opposite extreme cases that are beyond this range. At one extreme, 
in the most strongly "infra red" divergent hyper- membrane case p = q — 1, the formula ( |111| ) 
needs a modified factor (p — l)/2p. However, as discussed above, the quantities involved are, 
in this case, too highly cut off dependent for the result in question to have much significance, 
and anyway, as remarked above, their ultraviolet misbehaviour is so mild that hyper membranes 
can be satisfactorily treated without recourse to regularisation, so no analogue of the gradient 
formula is needed. 

On the other hand, at the opposite extreme, in the most strongly "ultra violet" divergent 
case, namely that of a simple point particle with p = 0, an appropriate regularised gradient 
formula is indeed something that is obviously needed. However, for a static configuration, the 
possibility of curvature does not arise at all in this "zero-brane" case, so that - while there is 
no reason to doubt the conjecture that it should still apply - there is no short cut whereby the 



formula (111) can be derived without further work, which will require recourse to the technical 
complications (due to the dimension sensitive nature of the time-dependent Green functions |T^ 1) 
involved in a fully dynamical treatment. 

Appendix: Integral formulae and notation 

The standard definition of the Euler Gamma function is provided by the integral formula 

T{z}= / e~H z ~ l dt. (112) 
Jo 

When the argument has an integer value n it can be shown, using the obvious recursion relation 
r{n + 1} = nT{n}, that it satisfies the well known relations 

r> + l} = n!, r{n + i} = M|0F. (113) 
In this article, we frequently need to evaluate integrals of the form 

'*W„ jiT?j* if ' (114) 

Whenever p > 1 and q > p, this integral will be convergent and will be expressible - via a change 
of variables t = 1/(1 + y 2 ) - in the form 

where B{a, b} is an Euler integral of the first kind - also known as the Beta function - that is 
specified by the formulae 

B{a, b} f t»-\l - tf-^t = ^§ . (116) 
Jo T{a + b) 

The values, as given by (|l~2|), of the surface area fi' 9 11 of a unit g-sphere are related for different 
values of q by the formulae 



n lp - 1] n lq - p - 1] " r{|} 1 2' 2 



B\V-^\, (H7) 
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and thus satisfy the useful relation 

q (p q-p } tt lp ~ 1] _ q-p 

2 B \2— + l ]W^-W^- (118) 



In order to facilitate the reading of the article, we conclude by providing the following table 
summarising the the notation used for the profile, Green, and other related functions in the three 
cases considered in this article. 





General notation 


Infinitely thin limit 


canonical case 


Potential 




<t> 




Source term 




P 




Profile distribution 


v u 


5 [q] 


8 


Green function 


-L— [<2,J>] 


G w 


% 

r W>J>] 


Radial profile function 




Radial Green function 
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